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ABSTRACT 

Context. The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key 
initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have 
focused on single stars. Most O stars are however found in multiple systems. 

Aims. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distribu¬ 
tions of binary sub-populations with one another as well as with that of presumed single stars in the same region, we aim to constrain 
the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the 
massive stars spin rates. 

Methods. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) 
to establish the projected equatorial rotational velocities ( v e sin i) for components of 114 spectroscopic binaries in 30 Doradus. The 
Vq sin i values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. v e sin i calibra¬ 
tion that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results. The overall v e sin / distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low- 
velocity peak (at sin i < 200 km s -1 ) and a shoulder at intermediate velocities (200 < i> e sin i < 300 km s -1 ). The distributions of 
binaries and single stars however differ in two ways. First, the main peak at v Q sin i ~100 km s -1 is broader and slightly shifted toward 
higher spin rates in the binary distribution compared to that of the presumed-single stars. This shift is mostly due to short-period 
binaries (P orb < 10 d). Second, the v e sin i distribution of primaries lacks a significant population of stars spinning faster than 300 
km s -1 while such a population is clearly present in the single star sample. The v e sin i distribution of binaries with amplitudes of 
radial velocity variation in the range of 20 to 200 km s -1 (mostly binaries with P orb ~ 10 - 1000 d and/or with q < 0.5) is similar to 
that of single O stars below v e sin i < 170 km s -1 . 

Conclusions. Our results are compatible with the assumption that binary components formed with the same spin distribution as single 
stars and that this distribution contains few or no fast spinning stars. The larger average spin rate of stars in short-period binaries may 
either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars 
and long period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of 
SB2 systems with P orb < 10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries of 
which the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most 
stars with i> e sin i > 300 km s -1 in the single star sample are actually spun-up post-binary interaction products. Finally, the overall sim¬ 
ilarities (low-velocity peak and intermediate velocity shoulder) of the spin distribution of binary and single stars argue for a massive 
star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a 
binary system. 

Key words, stars: rotation - stars: binaries: spectroscopic - Magellanic Clouds - Galaxies: star clusters: individual: 30 Doradus - 
line: profiles 


1. Introduction 

The spin rate that a massive star acquires at birth carries infor¬ 
mation about the poorly constrained processes involved in its 


* Based on observations collected at the European Southern 
Observatory under program ID 182.D-0222. 


formation (e.g., [Zmnecker &~^Y orke 2007] |Rosen et~aLl|2012| 


Krumholz 2014|. Once assembled, the spin of a massive star in 


fluences its nucleosynthesis and lifetime and more generally its 


evolutionary scenario (e.g 

Maeder & Meynet|2000; Brott et al. 

2011; Ekstrom et al. 2012 

Georgy et al. 2012) and supernova 

characteristics (e.g., Yoon & Langer 2005; Woosley & Heger 
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|2006HYusof et al.|2013[|Groh et al.|2014] ). For a review of these 
topics we refer to [Danger ([2012 ). 

Up to now, studies of the spin properties of massive stars 
have focused on the presumed-single star population with only 
some report i ng measurements of spectroscopic binaries (e.g., 
Penny 1996; Penny & Gies 2009). Here we present an analysis 


Agudelo etal. ( |2013| ) found a v e sin i distribution rotational peak¬ 
ing at about 80 km s _1 . After deconvolution assuming arbitrary 
orientations of stellar spin vectors this peak shifts to about 100 
km s“\ corresponding to -20% of the break-up speed. Because 
spin-down throu gh ste llar winds is limited for the vast majority 
of the sample, [Ramirez-Agudelo et ah] ( 201 3| ) suggested that this 
peak was representative of the initial distribution of spin rates of 
single stars, i.e., of the outcome of the formation process. 

[Dufton et al.| ( [2013| investigated the VFTS spin distribution 
of non-supergiant early-B stars and found a strong bi-modality. 
The distribution presents a low-velocity peak at < 100 km s -1 
and a broad and even stronger high-velocity component extend¬ 
ing from v Q « 170 to 350 kms -1 . The low- and high-rotation 
velocity regimes in the B-star sample are separated by a clear 
v G sin i gap extending from 100 to about 170 km s -1 . The strong 
bi-modality of the spin-rate distribution seen in the VFTS B-star 
sample is generall y not observed in the VFTS single O-type stars 
( [Ramirez-Agudelo et al.|2013| . A hint of bi-modality may how¬ 
ever be present in the 09 distribution and may suggest a smooth 
transition from the O star to the B star regimes. The weakly 
pronounced second O-star peak is however located at a slightly 
lower projected spin velocity and is not as broad as the one seen 
in the B-star distribution. Further data would also be needed to 
validate the statistical significance of such a peak in the 09 star 
population. 

According to the authors, the low spin velocity peak in the B- 
star distribution may be the result of magnetic spin-down, with 
long-lived magnetic fields possibly been generated through stel¬ 
lar mergers in close binary systems. Since close binaries are un¬ 
likely to contain stellar merger products, such a scenario would 
predict that the fraction of binary components with the lowest 
spin rates is reduced in the spin distribution of binary compo¬ 
nents. 

The observed v Q sin i distribution of O stars further presents 
a declining high-velocity tail at v e sin i > 300 km s -1 . By com- 
parison with the results of recent population synthesis ( |de M ink 


|et al.||2013| , [Ramirez-Agudelo et al.| ( |2013| ) suggested that the 


of the projected equatorial rotational velocity (v e sin i) properties 
of the O-type spectroscopic binary population observed in the 
context of the VLT-FLAMES Tarantula Survey (VFTS; Evans] 
et al.|2011b| ), a multi-epoch optical spectroscopic campaign tar¬ 


geting 360 O-type objects and hundreds of cooler stars in the 
30 Doradus region in the Earge Magellanic Cloud. Overviews of 
the VFTS project and the results obtained so far are provided in 
e.g.,|Evans et al.|(|2011a|),|de Koter et al.|(|2011||2013|) and|Sana| 
J20I3E 

For the single O -type stars in the VFTS survey, |Ramirez-| 


The O-type binaries in the VFTS are identified by their RV 
variations ( [Sana et ah 2013) and hence will be dominated by 
pre-interaction systems ( de Mink et al.|2014] ). To first order, this 
sample provides thus a view on the rotation properties of mas¬ 
sive stars which is unaffected by strong binary interactions, ei¬ 
ther through mass transfer or merging. In this respects, results 
from the binary sample may thus be confronted to that from 
presumed-single stars, where binary interaction effects are ex¬ 
pected to lead to strong signals (de Mink et al.|2013||2014| ). 

Furthermore, we may look at the role of tidal interactions 
in systems with periods on the order of 1-10 days. This tends 
to synchronize the rotation of the components with the orbital 
motion, leading to a spin-up of stars that at first have modest ro¬ 
tation rates, and to the synchronization of the angular velocities 
of primaries and secondaries in the same binary (jde Mink et al. 
12013). 

Finally, we may take out the strong binary effects from the 
single star spin distribution, that of tides from the binary compo¬ 
nent spin distribution, and then relate the differences, if any, in 
the two resulting distributions to differences in the star formation 
processes of single and binary stars. 

This paper is organized as follows. Section [2] introduces our 
sample. Section [3] describes the adopted methodology to esti¬ 
mate the projected spin rates of the components of the VFTS 
O-type spectroscopic binaries. Section]?] presents and analyses 
the obtained v e sini distributions. Our results are discussed in 
Section [5] and our conclusions are summarized in Section [6] 

2. VFTS data and O-type binary sample 


high rotational velocity tail may be dominated by spun-up post¬ 
interaction objects. Such binary evolution products are indeed 
expected to be seen as single stars in radial velocity (RV) sur¬ 
veys due to the absence of large RV variations, either because 
of coalescence of the binary components or because the mass 
donor has lost so much mass that the induced RV variations of 
the mass gainer - that dominates the system optical brightness 
after mass transfer - are below any practical detection limit for 
spectroscopic binaries ( [de Mink et al.||20l4| ). The mass gainer 
may also be truly single in cases where the binary is disrupted 
by the supernova explosion of the primary component. 


The VFTS campaign has been described in E vans et al. ( 2011b ). 
All spectra discussed here were obtained with the Medusa fibre- 
feed to the Giraffe spectrograph. The VFTS Medusa sample con¬ 
tains 332 O-type objects. Spectral classifi cation and RV mea- 
surements were presented in |Walborn et al. ( 2014} ) and Sana et al. 
( 2013| ), respectively. 

The method adopted for RV measurements was based on 
Gaussian fitting and, for each O star in the sample, pro¬ 
vided Doppler shift, amplitude and full width at half maximum 
(FWHM) values for a subset of Hei and Hen lines. 116 objects 
were classified as spectroscopic binaries (SB) on the basis of 
their statistically significant and large-amplitude RV variations 
(ARV >20kms _1 ). 

The identified O-type binary sample contains 85 single-lined 
(SB1) and 31 double-lined (SB2) spectroscopic binary systems. 
Figure [I] shows the spatial distribution of the binary stars. The 
systems are concentrated in the two clusters, NGC2070 and 
NGC2060 (6.7’ south-west of NGC2070), although a sizeable 
fraction is distributed throughout the field of view. 

3. Projected rotational velocities 

The projected rotational velocity, v e sin /, of stars can be mea- 
sured fro m the broadening of their spectral lines ( Carroll|1933[ 
Gray 1976] ). For OB-type stars commonly used methods include 


the direct measurement of the full width at half maximum (e.g., 
|Slettebak et al.|1975||Herrero et al.|1992[|Abt et al.|2002| ), cross¬ 
correlation of_theobserved_spectrum against a template spec¬ 
trum (e.g., Penny 1996 Howarth et al. 1997| ), and comparison 


with synthetic lines ca 


culated from model atmospheres (e.g., 


Mok iem et al.|2006| ). |Sim6n-Diaz & Herrero| ( |2007| 2014| ) apply 


a Fourier transform (FT) and a line-profile fitting method. 

The phase coverage of the VFTS binaries does not allow us 
to disentangle the spectra of the binary components, hence one 
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Fig. 2. v e sin i derived by means of Fourier transform or line-profile fitting methods (vertical axis; see Ramirez-Agudelo et al. 


versus FWHM for Hei /14387, 4713, 4922, and Hen 44541 (horizontal axis). The lines show the fitted FWHM - v e 
relations described by Eq.[l]and Table [T] 
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Fig. 1 . Spatial distribution of the O-type binary systems in the 
30 Dor field of view. Dots indicate SB1 systems (85 objects) 
and crosses SB2 systems (31 objects). The dashed circles have 
radii of 2.4' and are centered on the central cluster NGC 2070 
and on NGC 2060, at 6.7’ to the south-west of NGC 2070. (A 
color version is available in the online version of the paper). 


cannot easily apply FT or line profile fitting methods. To ob¬ 
tain the projected rotational velocities of the binary sample in 
a homogeneous way, we thus focus on the FWHM of specific 
spectral lines. The FWHMs are indeed easily quantifiable and 
less affected by line blending and/or the presence of faint or 
undetected companions than more sophisticated methods. The 
drawbacks are relatively large error bars on individual mea¬ 
surements (see below) and a limited discriminative power at 
v e sin i < 100 km s -1 due to the fact that the method does not 


disentangle rotational broadening from other types of broaden¬ 
ing such as macro-turbulence and Stark broadening. 


3.1. FWHM - v c sin i calibration 

For the single O-type stars in the VFTS sample, we used the pro¬ 
jected rotational velocities previously obtained using both a FT 
and a line profile fitting method ( [Ramirez-Agudelo et al.|2013| . 
The FT method is able to disentangle rotational broadening from 
other broadening agents (such as for instance broadening due to 
macro-turbulent velocity fields) in cases where rotational effects 
dominate. For stars where rotation does not dominate the broad¬ 
ening, the line fitting method also takes into consideration the 
presence of additional broadening (for more details see Ramirez- 
| Agudelo et al.|2013|). 

The Hei 44387, 4922, 4713, and Hen 44541 lines studied 
Ramirez-Agudelo et al.|(|2013) have also been used for RV 


measurements ( [Sana et al.|2013| and thus have both v e sin i and 
FWHM measurements available. This allows for a calibration of 
the FWHM values in terms of v Q sin i. 

Figure [2] compares v Q sin i vs. FWHM for the helium lines 
listed above. The number of stars for which this comparison can 
be made differs per spectral line and depends on spectral type 
and data quality. The main reason for this is that the lines need 
to pass a set of quality criteria to ensure reliable measurements, 
i.e., they should be of sufficient strength relative to the signal- 
to-noise ratio of the spectrum and to have a sin i measurement 
that is above the resolution limit ( [Ramirez-Agudelo et al. |20 1 3) . 

All four lines show a clear trend between v e sin i and FWHM. 
The scatter in the relation for He n is considerably larger at rela¬ 
tively low v Q sin i than for the three Hei lines. This results from 
the fact that the Stark broadening is more pronounced in He n 
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Table 1. FWHM - v c sin i calibration based on the VFTS single 
O-type star sample. The columns list the diagnostic line (column 
1), the number of stars for which the line fulfills the selection 
criteria (column 2), the coefficients of fit function Eq.[l](columns 
3 and 4) and the root mean square (rms) deviation between fit 
and data in km s -1 (column 5). 


Diagnostic line 

# stars 

a 

b 

rms 



(km s -1 /A.) 


(km s -1 ) 

He 1 44387 

102 

28.99 + 1.55 

1.30 + 0.03 

25 

Hei 44713 

60 

35.51 + 1.81 

1.20 + 0.04 

15 

Hei44922 

69 

26.00 + 2.01 

1.26 + 0.04 

25 

Hen 44541 

92 

10.17 + 1.27 

1.69 + 0.07 

35 


lines , and that a_treatment of_this broadening is not properly 
included in |Ramfrez-Agudelo et al. ( 2013| . 

To calibrate FWHM as a function of v Q sin i we obtain a phe¬ 
nomenological description of the trends seen in Fig. [2] We suc¬ 
cessively adjusted polynomials of order 1, 2 and 3 as well as a 
power-law function with and without an offset. We perform F- 
tests to decide whether additional fit parameters result in signifi¬ 
cantly better representations of the data. All fit functions, except 
the first order polynomial, give similar quality results. As the 
power law without an offset has the fewest free parameters, we 
opted for this representation i.e., 


v e sin i = ax FWHM , 


( 1 ) 


where a and b are the fitting parameters. 

Table Q] summarizes the best fit values and their uncertainties 
for each of the lines. The root mean square (rms) residuals are 
about 15 to 25 km s -1 for the Hei lines and reach 35 km s -1 for 
the He n 44541 line. In the next sub-section, we lay out our strat¬ 
egy in using these four relations and investigate the uncertainty 
in the v e sin i determinations that follow from this approach. 


3.2. FWHM- based v t sin i values for the single star sample 


We first apply the FWHM-t> e sin i calibration to the sample of 
216 presumably single O-type stars, as this provides an internal 
consistency test. For four stars (VFTS 400, 444, 565, and 587), 
no FWHM measurements of any of the four diagnostic lines are 
available. These have been excluded from our comparison, leav¬ 
ing 21 2 stars. 

In Rarmrez-Agudelo et al. ( |2013| ) the diagnostic lines have 
been divided into two groups, one of relatively good quality 
lines (Group A) and one of lesser quality lines (Group B). The 
Group A lines include the Hei lines at 4387, 4713, and 4922 A. 
Only if these lines are very weak or not present, v e sin i relied on 
Hen 44541. These latter measurements are referred to as Group 
B. We also adopt this division of measurements in two groups 
in order to decide which lines provide the most reliable value 
of v Q sin i. If more than one He i line is available in a Group A 
source, we use the unweighted average to obtain the mean pro¬ 
jected rotational velocity. As this mean is based on at most three 
values, the standard deviation tends to underestimate the true dis¬ 
persion. We thus adopt the rms residuals of Table [I] as error bars 
for our measurements, unless the standard deviation of a given 
star is larger. In that case we adopt the computed value. 


1 While an extra offset would formally express the fact that spectral 
lines have an intrinsic width - hence would be more physical - adding 
this extra parameter yields no improvements on the fit quality and we 
opted for the simpler function. 



Fig. 3. Comparison of v e sin i values derived from the FWHM - 
Vq sin i calibration (Eq. [TJ and those of Ramirez-Agudelo et al. 
( |2013| ), which are based on more sophisticated methods. The 
dashed line shows the one-to-one relation. 


Figure 3] compares the v Q sin i values based on the FWHM 


v e sin i calibration with those obtained in Ramirez-Agudelo et al 


(2013). The systematic difference is less than 5 km s -i and the 
standard deviation is about 40 km s -1 . For slowly spinning stars, 
the calibration presented here tends to somewhat overestimate 
v e sin i. Indeed in this low-velocity regime the FWHM-t> e sin i 
calibration fails to disentangle the contribution of rotat iona l 
broadening from other broadening mechanisms (see Sect. 

At large velocity (t> e sin i > 400 km s ), our FWHM-based mea¬ 
surements tend to underestimate v e sin i compared to values ob¬ 
tained with the FT or profile fitting method. It may suggest the 
presence of a higher-degree component in the calibrations but 
our previous tests show that, if present, the latter could not be 
firmly identified considering our data quality. 

Figure [4] compares the v t sin i probability density distribu- 
tio n (top pane l) and v Q sin i frequency distribution (bottom panel) 
of [Ramirez-Agudelo et al.| ( |2013| ) with that derived from the 
FWHM- Vq sin / calibration, ad opting a bin size of 40 kms -1 
as in [Ramirez-Agudelo et al.| ( [2013 ). The first bin of the new 
calibration is depopulated relative to the Ramirez-Agudelo et al. 
( 2013[ ) analysis because of the neglect of additional broadening 
processes. Such extra broadening may affect v Q sin i for stars ro¬ 
tating below 150 km s -1 , and also explains the shift of the peak 
of the distribution to somewhat higher v e sin i values. At veloc¬ 
ities above 200 km s -1 , the distribution is not affected in a sta¬ 
tistically significant way, and the small (one bin) changes that 
occur can be explained by the random uncertainties in the two 
methods. 

Importantly, the general shape of the distribution - a peak at 
about 100 km s -1 and a high-velocity tail - is preserved using 
the FWHM-f e sin i method. Despite relatively large individual 
measurement errors, this method still allows us to statistically 
investigate the rotational properties of massive stars, including 
those of binaries for which more sophisticated methods are more 
cumbersome to apply. 


4. The Vq sin i distribution of spectroscopic binaries 

In this section, we construct and investigate the v & sin i distri¬ 
bution of stars in O-type binary systems. Out of 116 primaries 
Vq sin i values could not be established in two cases, VFTS 531 
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Projected rotational velocity (km/s) 



Fig. 4. Cumulative (upper panel) and frequency (lower panel, 
with Poisson error bars) distributions of the projected rotational 
velocities of the presumed-single O-type stars obtained usin g 
the FWHM-t> e sin i calibration (Eq. [TJ and those of Ramirez- 
|Agudelo et al. ( |2013| ). 


and VFTS 810, as th e spectra do not pass the quality criteria de¬ 
scribed in Sect. 13. H and thus these two stars were excluded from 
our analysis. 

From here on, we only use v Q sin i values obtained using the 
FWHM-f e sin i calibration, i.e., also for the set of single stars. 
This allows for a direct and consistent comparison of the single 
and binary samples. 

We further perform statistical comparisons of the v Q sin i dis¬ 
tribution from different sub-samples in order to identify statis¬ 
tically significant signals in our data. We use Kuiper’s test of 
hypothesis (KP, |Kuiper|p~960| to search for significant differ¬ 
ences between the parent populations of the considered sam¬ 
ples. Results are reported in Table [2j where we list the sum D of 
the absolute values of the most positive and most negative dif¬ 
ferences between the two cumulative distribution functions that 
are being compared and the Kuiper probability /?k- Specifically, 
KP tests allow us to test the null hypothesis that two observed 
distributions are randomly drawn from the same parent popula¬ 
tion. This test is similar to the more well-known Kolmogorov- 
Smirnov test but presents the advantage that it is as sensitive to 
variations in the tail as it is to variations in the middle of the 
distributions. 

In Table [3] we also present the frequency of stars in a given 
sub-population that are spinning faster than the specified limit. 


Table 2. Kuiper’s test statistics. The figure index is given in 
the first column. The considered samples and number of stars 
in each of these samples are listed in Cols. 2 to 5. Column 6 
specifies the KP test statistic D and column 7 the probability p k. 


Fig. 

Sample 1 

Sample 2 

nl 

n2 

D 

Pk 

(%) 

5 


Presumed-single 

Primaries 

212 

114 

0.20 

3 

7 


Primaries 

Secondaries 

114 

31 

0.18 

95 

9 


Low-ARY 

High-ARY 

85 

29 

0.24 

57 

8 


Primaries (SB2) 

Secondaries (SB2) 

31 

31 

0.16 

99 



For Vq sin i < 170 km s 

i 




10 

Presumed-single 

Low-ARY 

156 

64 

0.15 

76 


Table 3. Frequency of stars (/) from different sub-samples that 
are spinning more rapidly than the specified limit. 68% confi¬ 
dence intervals (approx. ± 1 sigma error bars) were computed 
following binomial statistics ( Sana et al. |2014| ). 




/Oe sin 0 


Sample 

> 200 km s 1 

> 250 km s -1 

> 300 km s 1 


Presumed-single 

0.193 

± 

0.028 

0.156 

± 

0.024 

0.104 

±0.019 

Primaries 

0.175 

± 

0.035 

0.088 

± 

0.026 

0.026 

±0.018 

Primaries (SB2) 

0.161 

± 

0.065 

0.032 

± 

0.032 

0.000 

(< 0.032) 

Secondaries 

0.097 

± 

0.065 

0.032 

± 

0.032 

0.032 

± 0.032 

Low-ARV 

0.141 

± 

0.035 

0.082 

± 

0.035 

0.035 

± 0.024 

High-ARY 

0.276 

± 

0.069 

0.103 

± 

0.069 

0.000 

(< 0.034) 


We will use these statistics to assess the significance of the pres¬ 
ence of a high-velocity tail in these sub-populations. 

4.1. Primaries vs. single stars 

4.1.1. Observed distributions 

Figure [5] compares the v e sin i distributions of presumed-single 
stars and primaries. There is a surprising overall qualitative 
similarity between the two distributions. Both are dominated 
by a low-velocity peak (< 200 kms -1 ) which contains ~80% 
of the samples, and show a shoulder at intermediate veloci¬ 
ties (200 kms -1 < L> e sin/ < 300 kms -1 ). The distribution of 
single stars also features a clearly identifiable high-velocity (> 
300 km s -1 ) tail, which is hardly seen in the distribution of pri¬ 
maries. 

Though qualitatively similar the main low-velocity peaks are 
not identical. The main peak in the distribution of the primaries 
is wider than that of the presumed-single sample and overpopu¬ 
lates the distribution in the region between 100 and 300 km s -1 
by 18%. At large projected spin rates (t> e sin/ > 300 kms -1 ), 
there is a deficiency of very fast rotating primaries with respect 
to the presumed-single stars. The presumed-single sample har¬ 
bours 22 stars out of 212 that exhibit projected rotational veloc¬ 
ities larger than 300 km s -1 (10.4 ± 1.9 % of that sample) while 
the sample of primaries only has 3 stars out of 114 (2.6 ± 1.8 %, 
see Table [3]). A KP test confirms that parent v c sin i distribu¬ 
tions of primaries and of single stars are statistically different 
(Pk = 3%). 


4.1.2. Impact of inclination distributions 

Binary systems that are seen closer to edge-on, i.e., that have 
a large inclination angle, are more easily identified as binaries 
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Table 4. Baysian estimates of the rotation distribution parame¬ 
ters of Eq.[2] Reported values correspond to the 50th percentiles 
and confidence interval to the 16th and 84th percentiles of the 
posterior distributions. 


Parameters 

Primaries 

Single stars 

Adjusted parameters 

a 

7.83-JH 

8.94!°;™ 

p 

0-0476" 

n n79Q +0 - 0086 
u.u / _ 00121 

n 

215!!?? 

320!j° 

cr 

1423J * 1 

113+g 

Iy 

0-94!^ 

0.7its? 


Computed parameters 

H y = a 1/3 164+!! 

122+s 

°v = V«/^ 2 59 ~l 

41!? 


because of larger projected radial velocity variations. This may 
introduce an observational bias in the spin rate distribution of 
binaries relative to that of single stars if the spin axis of the com¬ 
ponents is preferentially oriented along the orbital rotation axis. 
For identical intrinsic spin distributions, the observed distribu¬ 
tion of single stars might thus be shifted to lower v e sin i than 
that of binaries, potentially explaining part of the shift that we 
see in Fig. [5] 

An estimate of the most extreme effect of this bias is pro¬ 
vided by comparing the deconvolved v e distribution of single 
stars (fig. 16 of Ramirez-Agudelo et aL||2013 ) to the observed 
v e sin i distribution of binaries in Fig. [5] The latter indeed corre¬ 
sponds to the intrinsic distribution of binaries under maximum 
detection likelihood corresponding to i = 90° for all systems. 
Even in these conditions, the main peak of the binary distribu¬ 
tion is still broader and shifted towards higher v Q compared to the 
main peak of the single star distribution. Furthermore, we show 
in Appendix [A] that the actual distribution of inclination angle 
of the binaries in our sample is not sufficiently different to that 
expected from random orientation in 3D space to significantly 
contribute to differences seen in Fig. [5] We conclude that the 
possible bias induced by a preferred binary detection for higher 
inclination system cannot be invoked to explain the difference 
of the v Q distributions in the primaries and presumed-single star 
samples. 


4.1.3. Intrinsic v e distributions 

In Appendix [A] we use a simulated distribution of orbital incli¬ 
nation of the binaries in our sample and a Bayesian method to 
compute the intrinsic distribution of rotational velocities of the 
primaries in our sample. We adopt a general analytical form of 
the rotation velocity probability density function /(t> e ), s imilar to 
the one found for single O stars within VFTS ( [Ramirez -Agudelo' 
et al.|2013| , i.e., a combination of a gamma distribution and of a 

normal distribution: 


f(Pe) = ly g(v e a,(3) + (1 - Iy) N(v e ‘,fl, O 2 ) 


where 
g(x; a,/3) 


N(x; n, cr 2 ) 


pa 

Off) 




1 + erf(/j / V2<x) V2 ncr 


( 2 ) 

(3) 

(4) 




Projected rotational velocity (km/s) 


Fig. 5. Cumulative (upper panel) and frequency (lower panel, 
with Poisson error bars) distributions of the projected rotational 
velocities of the presumed-single O-type stars (y e sin i - FWHM 
method) and O-type primaries. (ARV < 20 km s -1 ) 


In Eq. [4j erf is the error function. It is used to renormalize the 
Gaussian component to unity between x = 0 and oo. Similarly, I y 
and (1 - Iy) are the relative contributions of both distributions to 
/(v e ), such that /(v e ) is also normalized to unity when integrated 
over positive v e . Such renormalizations were not implemented in 
the same way in Ramirez-Agudelo et al. j2013| . For consistency 
in the comparison, we thus also applied the Bayesian analysis 
method to the single star y e sin i measurements obtained through 
our FWHM calibration curves. 

Table [4] provides the 50th percentiles as well as the 68% 
confidence interval of the marginalized posterior distributions 
of each of the parameters of Eq. [2] It also provides the mean 
(piy = a IP) and standard deviation (<x r = yjajfp) of the gamma- 
component in the rotation distributions. Projections of the cor¬ 
responding ma rginal ized post erior probability distributions are 
shown in Figs. |A.3| and A.4 Finally, Fig. [6] compares an his¬ 
togram of the observed v e sin i data with the obtained best-fit ro¬ 
tational velocity distributions (see Appendix |A| for details). 

For the single star distribution, the best-fitjparameters are 
different from those obtained by [Ramirez-Agudelo et al. ( 2013| ). 
This likely originates from the different analytical forms and fit¬ 
ting methodology used. It may also have origins in the derivation 
of the Vq sin i measurements (now based on FWHM for consis¬ 
tency). Still, the overall shape of the recovered intrinsic rota¬ 
tional velocity distribution is in good agreement (see Fig. |A.2| ), 
which allows us to validate the method. 
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The Bayesian analysis allows us to confirm and quantify our 
previous conclusions: 

(i) The peak of the intrinsic rotation distribution (y e ) is shifted 
by almost 50 km s -1 towards larger values in the sample of 
primaries; 

(ii) The main peak is 50% broader in the distribution of pri¬ 
maries (cr y = 59 km s -1 compared to 41 km s -1 in the single 
star distribution); 

(iii) The high-i; e sin i tail is very limited or altogether absent in 
the sample of primarie^] 

Our analysis is the most complete, bias-corrected and con¬ 
sistent comparison of the rotation properties of single and bi¬ 
nary massive stars. It reveals small but fundamental difference 
between both samples, which we will discuss in Sect. [5] We also 
derive a suitable analytical representation of the intrinsic dis¬ 
tribution of rotation rates that can be used, e.g., in population 
synthesis. In doing so, two aspects have to be kept in mind: 
(i) the derived distributions are present-day distributions, not 
initial ones; (ii) the lack of extremely slow rotators (t> e sin/ < 
40km s -1 ) in the distributions results from the limitation of our 
Vq sin i measurement method. Higher resolution data and, for the 
binaries, spectral disentangling, would be needed to further our 
analysis at low spin rates. 

4.2. Secondaries vs. primaries 

In Fig. [7] we show the distributions for primaries (114 stars) and 
secondaries (31 stars). Qualitatively, there is no clear difference 
between both distributions, i.e., the peak is at a similar rotational 
velocity and both distributions are equally broad. The primary 
distribution presents two stars with v Q sin i > 400 km s -1 . These 
two outliers represent less than 2+2 % of the total primary sam¬ 
ple and do not impact a comparison of the two distributions. A 
KP test shows that both distributions are compatible with one 
another ( p K = 95 %). 

Similarly, Fig. [8] displays the probability density distribu¬ 
tions of the primaries and of the secondaries in our sample of 
31 SB2 binaries. These systems mostly have large radial veloc¬ 
ity variations (ARV > 200 km s -1 ) indicating that the population 
is dominated by close binaries. A Kuiper test yields D = 0.16 
and pk = 99%. This is thus consistent with the hypothesis of 
tidal synchronization, assuming the sizes of the accompanying 
primaries and secondaries are comparable. 

4.3. Small vs. large RV variation systems 

The size of the sample of primaries allows us to explore y e sin i 
distributions for sub-populations. By selecting a limit for the am¬ 
plitude of the radial velocity variations ARV ^ t of 200 km s -1 
we divide the primary sample into two sub-samples, i.e., stars 
that display ARV < ARVi imit (henceforth low-ARV) and stars 

2 The Gaussian-component - used to add enough degrees of freedom 
in the fit to be able to model a possible high-u e sin i tail - is not sig¬ 
nificant, at the 5% significance level, in the distribution of primaries. 
This yields very broad marginalized posterior distributions for p and 
cr, hence poorly defined parameters for the Gaussian component as 
indicated in Table [4] Within the 68% confidence interval, the ampli¬ 
tude of the Gaussian component does however reach 1 - l y - 0.21, in 
which case one could consider it to be important. However, inspection 
of Fig. | A.3 1 reveals that the parameters of a Gaussian that would have 
1 -ly « 0.21 are given by p « 200 and cr & 50, i.e., such a Gaussian does 
not contribute significantly to a high-i> e sin i tail above « 250 km s -1 . 



(Projected) rotational velocity [km s 1 ] 



(Projected) rotational velocity [km s 1 ] 


Fig. 6. Rotational velocity distributions for 50 random samples 
of the MCMC chain for the primaries (upper panel) and the sin¬ 
gle stars (lower panel). The intrinsic rotational velocity distribu¬ 
tions are shown in red. In blue we show the corresponding distri¬ 
butions projected (assuming the sin i distribution inferred for the 
detected binaries in VFTS for the primaries and a random ori¬ 
entation of the system in 3D space for the single stars) and con¬ 
volved with a Gaussian error distribution with cr v = 20 km s -1 . 
Although the observed data is represented here by a histogram, 
the fit to the data is performed discretely, i.e., without any bin¬ 
ning. The histogram bins are chosen according to Knuth’s rule 
( |Knuthl2006| ). 


with ARV > ARViimit (high-ARV). The latter are systems for 
which tidal synchronization is expected to be important during 
the main-sequence phase. Indeed, given our sample prop erties 
and using Monte Carlo simulations described in Appendix | A. 2| 
we can show that 88% of systems with ARVii m i t > 2 00 km s~ J 
have orbital periods less than 10 days (see section [571]) . Most of 
the systems in the low-ARV sub-sample will not suffer from tides 
-75% of them are expected to have P > 10 d and/or q < 0.5. 
The exception is of course for short-period similar-mass systems 
that are seen at low inclination angles or for which the RV-curve 
has been poorly sampled. 

Figure [9] plots the v e sin i distributions of the low-ARV (85 
stars) and high-ARV (29 stars) sub-samples. At low rotation ve¬ 
locities the distribution of high-ARV, compared to the low-ARV, 
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Fig. 7. Cumulative (upper panel) and frequency (lower panel, 
with Poisson error bars) distributions of the projected rotational 
velocities of the O-type primaries and secondaries. 


Fig. 8. Probability density distribution of u e sin i for the 31 SB2 
systems, for which v e sin i of both the primaries and secondary 
component can be constrained. 


appears shifted by one bin (« 40 kms -1 ) to higher rotational 
velocities. The high-velocity tail (i.e., v e sin i > 300 km s -1 ) of 
the low-ARV distribution has three primaries (3.5+2.4 % of that 
sample), while there is no star in the high-ARV distribution. 
Despite these apparent differences, a KP test does not allow us 
to reject the hypothesis that both samples are drawn from the 
same parent distribution ( p K = 57 %), possibly as a result of 
the limited size of the high-ARV sample. The weighted mean of 
the samples (130+2 km s -1 and 153+5 km s -1 for the low- and 
high-ARV samples, respectively) are however significantly dif¬ 
ferent from one another, confirming a higher average spin rate of 
the high-ARV (short periods) systems. 


low-velocity peak to be the representative of the outcome of the 
formation process where in (projected) rotational velocity it ex¬ 
tended up to ~ 160 km s -1 . Indeed, for v Q sin i < 170 km s -1 both 
distributions are very similar and a KP shows that both distribu¬ 
tions are compatible with each other ( p K = 76 %, see Table [2]). 

For the binaries with projected rotational velocities larger 
than 170 km s -1 both distributions behave differently. The sam¬ 
ple of presumed-single stars presents a significant contribution 
of stars rotating more rapidly than 300 km s -1 , while the wide¬ 
binary sample does not. 

5. Discussion 


4.4. Single stars versus /otv-ARV binaries 


We compare the sub-sample of low-ARV sources with that of the 
presumed-single stars in Fig. [TO] As pointed out in Sect. |4.3| the 
primary stars in the sample of long period systems are not likely 
to interact with their companion during their main-sequence 
lives (where we have used a low-ARV as a proxy for long pe¬ 
riods, barring low inclination systems and orbit sampling ef¬ 
fects; see Sect. |4.3| ). As the bulk of the binary systems and the 
presumed-single systems are on the main sequence, a compar¬ 
ison of the two may thus help to address the question whether 
or not the primaries in wide (enough) binaries and presumed- 
single stars leave the formation process with similar spin-rate 
distributions. In Ramirez-Agudelo et al.[ ( |2013 ) we claimed the 


The dominant feature in the u c sin i distributions of the O-type 
primaries and secondaries in 30 Dor is their low-velocity peak. 
Compared to the rotational distribution of single stars, there is 
a low occurrence of stars spinning faster than 300 km s -1 in the 
distribution of primaries. The well developed high-velocity tail 
seen in the presumed-single stars is thus nearly absent. In this 
section we consider possible physical mechanisms that may be 
impacting the global shape of the initial v e sin i distribution of 
binary components, notably tidal interaction and mass transfer. 


5.1. Signatures of tidal interaction 


As discussed earlier, the bulk of our binaryj>ample is expected 
to consist of pre-interaction systems, de Mink et al. (2013) study 
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Fig. 9. Cumulative (upper panel) and frequency (lower panel, 
with Poisson error bars) distributions of the projected rotational 
velocities of the O-type primaries with ARV < ARVnmit (low- 
ARV) and ARV > ARVi imit (high-ARV). 




Projected rotational velocity (km/s) 


Fig. 10. Cumulative (upper panel) and frequency (lower panel, 
with Poisson error bars) distributions of the projected rotational 
velocities of the presumed-single O-type stars and the sub¬ 
sample of low-ARV sources. 


the effects of binary evolution on the rotation rates of primaries 
and secondaries as a function of initial period and initial mass 
ratio. Prior to the primary star filling its Roche lobe, changes in 
stellar structure, mass loss and tides may affect the spin rates. 

Following Weidn er~& Vink| ( |2010| ), we estimate that about 
20% of the stars in our total sample have masses above 40 M 0 , 
with the most massive stars (accounting for ~4% of the total) 
being roughly 60-80 M©. It is these stars that develop winds 
that are strong enough to facilitate significant angular momen¬ 
tum loss during their main-sequence evolution (see Vink et ah] 


2001| ). As we focus on the ensemble properties, this implies that 


the bulk of the stars will not have stellar winds strong enough to 
spin them down significantly on the main-sequence timescale. 
Exceptions may be stars with strong large-scale magnetic fields, 
which we discuss separately below. 

Stellar expansion due to changes in internal structure will 
only have a minor effect on the spin rate when the star ap¬ 
proaches the end of its main-sequence lifetime, and can be fur¬ 
ther neglected ( |Brott et al.|2011||de Mink et al.|2013| ). The most 
important effect on the spin rate before Roche lobe overflow 
takes over is therefore that of tides. In systems where the radii 
of the comp onents is comparable to the separation between the 
stars < |de Mink et al.|2013| >, tidal interaction strives to synchro¬ 
nize the rotation of the binary components with the orbit. For 
massive binaries this corresponds to systems with orbital periods 
up to about a week. Once synchronized the rotational velocity of 


the co-rotating primary increases mainly because of secular stel¬ 
lar expansion. 

Like the single O-type stars in 30 Dor, the v c sin i distribu¬ 
tions of primaries and secondaries also present a low-velocity 
peak at around 100 km s -1 , though it is broader and skewed to 
somewhat higher spin rates (see Figs. Under the assump¬ 
tion that the initial spin distributions of binary components and 
single stars are the same, this difference is consistent with tidal 
interaction playing a role in our binary sample. Effects of tides 
and, in a minority of cases, possib le slow Case A mass transfer 
(see Figures 2 and 3; |de Mink et al.|2013| , appear most visible 
in the contrast between the projected spin-rate distributions of 
low-ARV versus high-ARV systems (see Fig. [9]), the latter corre¬ 
sponding to closer binaries and displaying larger spin rates (see 
Sect. @3}. 

To investigate the effect of tidal locking we make use of a 
diagram showing the relation between projected rotational ve¬ 
locity and amplitude of the RV variations (see Figs. ED an d[T2j>. 
The timescale of tidal synchronization is a function of the pri¬ 
mary mass Mi, the mass ratio of the primary and secondary 
q = M 2 1 and the orbital period R 0 rb- For our sample we lack 
information about the latter and therefore we use the amplitude 
of the radial velocity variations ARV as a proxy of the period. 

For a (hypothetical) fully synchronized binary system v e sin i 
and ARV are fully determined by R 0 rb> Mi, the radius R of the 
primary and q through Kepler’s law. The thick green line in 
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Fig. 11. Projected rotational velocity vs. amplitude of the radial 
velocity variations in the primary stars sample (114 stars). The 
shaded regions show approximately where we expect tidal syn¬ 
chronization to play a role during the main sequence evolution 
for a primary mass of 20 and 60 M© (with R = 8 and 10 R©, re¬ 
spectively). Here we varied the mass ratios q from 0.25 and 1 
and periods P 0 rb from 0.25 to 10 days. The labels on the top and 
right hand axes show the corresponding orbital and spin period 
for a typical circular binary system assuming a 20 Af© primary 
and a mass ratio of 0.5. 


Fig. |TT] indicates this relation for a typical system assuming a 
synchronized orbital and spin period between 10 days (lower 
left) down to the minimum period at which the primary star 
would be filling its Roche lobe (upper right). Here we adopted a 
primary mass of 20 A/©, with R = 8 R©, and a mass ratio q = 0.5. 
The shaded region shows the area spanned by similar relations 
varying the mass ratio between q = 0.25 (left side boundary of 
the green area) and q - 1 (right side boundary). The thick gray 
line and grey shaded area show the same adopting a primary 
mass of 60 Af© with R - 10 R Q . It stretches out to both higher 
rotation rates, because v e sin i scales with the stellar radius which 
is larger for more massive stars and to larger A RV which is pro¬ 
portional to Mj /3 . 

Tides are most effective for short periods systems ( |Zahn 
|1975[ |1977| with a weak preference for equal mass systems. 
Stars in the upper part of the shaded areas are expected to be 
fully synchronized as the timescale for synchronization becomes 
(much) less than 1 percent of the main-sequence lifetime ( fde] 
|Mink et ak]|2009[ ). For binary systems with periods larger than 
about 10 days and/or extreme mass ratios we do not expect tides 
to have a significant effect during the main sequence evolution. 

Forty VFTS primaries fall within the shaded zones. 
According to the above argument, they are expected to be syn¬ 
chronized or to tend toward synchronization. Given their v Q sin i 
properties, these systems are sufficiently numerous to explain the 
overpopulation of primaries in the 100-300 km s -1 range com¬ 
pared to the single star sample (Sect. |4.l| ). 

For the SB2 systems one can in principle check whether the 
two stars have synchronized. In Fig.[l2]we show the 31 stars for 
which we have obtained A RV and v Q sin i for both the primary 
and the secondary components. Most of the systems fall in the 
green zone. This is a selection effect as SB2s are most easily rec¬ 
ognized if the radial velocity shifts are large ( Sana et al. |201 1| ). 
All primary and secondary pairs have been connected with a 


Fig. 12. Projected rotational velocity vs. amplitude of the radial 
velocity variations for secondaries in SB2 systems (31 stars) to¬ 
gether with their corresponding primaries, linked with arrows. 
The definition of the green and gray regions and the labels on 


the top and right hand axes are as in Fig. 11 


straight line. If a pair is synchronized the angular spin velocities 
are the same. In that case one expects the lines to be horizontal in 
case of equal sized stars, and about horizontal when allowing for 
size differences. Indeed, by far the most systems show an about 
horizontal line connecting the two components and we thus in¬ 
terpret this property as the signature of tidal synchronization. 
Two distinct outliers are both located in the lower left corner of 
the diagram, where the timescale for tidal synchronization is the 
longest. 

Using population synthesis de Min k et ak| ( |2013| ) show that 
tides are expected to only lead to relatively few systems with 
velocities larger than 200 km s -1 , and that -300 km s -1 is the 
limit that tides can achieve. Close binary systems that potentially 
could have acquired spin rates much in excess of 300 km s -1 
at birth would be spun down to -200 km s -1 . Indeed, this ex¬ 
pectation too seems consistent with our findings: none of the 
components of the close systems (using ARV > 200 km s -1 as 
a proxy f or id entifying such systems) spin faster than this limit 
(see Sect. |5.3| ). Only when mass transfer occurs, later-on in the 
evolution, the secondary may reach spin rates of up to about 600 
km s -1 . 


5.2. Magnetic braking 

If a massive star possesses a sufficiently strong magnetic field, 
the interplay of the stellar winds with the stellar magnetosphere 
may lead to a rapid spin-down of the star. The efficiency of the 
spin-down depends both on the strength of the magnetic field 
and on that of the winds ( jUd-Doula et al.12009] ). For known mag¬ 
netic galactic O stars, the timescale of magnetic spin-do wn is of 
the order of several 0.1 Myr to 3 Myr (Petit et al . 2013] ). Given 
the lower metallicity, hence lower wind strengths of stars in the 
30 Doradus region, their spin-down timescales will be longer. 
Using Ud-Doula et al.| formalism, we estimated them to be of 
the order of 6 to 25 Myr for a 1 kG field, and a factor of ten 
shorter for a 10 kG field. 

Possibly up to 10% of the Galactic O stars host a long-lived 
large scale magnetic field ( [Morel et al.|2014[|~Wade et al.|2014 ). 
However, only one has a known magnetic field stronger than 
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10 kG (NGC 1624-2, [Wade et al.|2012||Petit et al.|2013> , a con- 
dition that seems required but rarely met for magnetic braking 
to efficiently affect LMC O stars. While the detected magnetic 
fields in Galactic OB stars are believed to be in a stable con¬ 
figuration, Fossati et al. (2015, in prep.) argue that the decay 
timescales in massive main-sequence stars may be shorter than 
their lifetimes. If the field decay timescale is longer than the 
magnetic breaking one (for stars to have time to spin-down), but 
fast enough for most strong magnetic fields to elude detection, a 
substantially larger fraction of stars may have been affected by 
magnetic breaking than what the current detection rate of strong 
magnetic fields suggests. As a consequence one needs to con¬ 
sider the possibility that magnetic braking may have some im¬ 
pact in our sample. 

The origin of the magnetic fields in massive stars is still un¬ 
der debate. One possibility is that of a fossil field that is frozen 
into a stable large-scale configuration during the star-formation 
process (Duez & Mathis 2010). Another is that the fields are 
formed during a strong binary interaction, e.g., a stellar merger 
(e.g., |Wickramasinghe et al.|2014 ). Merger products are not ex¬ 
pected to be found in close binaries, as then the pre-merger triple 
star would have been highly unstable. This agrees with the fact 
that the vast majority of magnetic massive stars seems to lack 
a close binary component, in contrast with their non-magnetic 
counterparts (Alecian et al. 2014| ). 

If magnetic spin-down is affecting our sample and if mag¬ 
netic fields are generated by binary interaction, one thus expects 
a reduced number of slowly spinning components in close bi¬ 
naries compared to single stars or to wide binaries. However, 
tidal interaction in close binaries leads to synchronization, such 
that this effect is not detectable. The resulting slowly spinning 
post-interaction products would mostly show up in the sample 
of presumed-single O stars. This would move the main peak 
in their initial spin distribution towards lower spin velocity. A 
similar effect may affect the wide binary distribution. Indeed, 
systems born as hierarchical triples may now appear as a wide 
binary, after the close pair has interacted. The post-interaction 
products may spin slowly because of magnetic spin-down. 

To conclude, this discussion illustrates the possibility that 
the observed shift of the spin distribution main peaks in the 
primaries and in the single stars samples may partly be due to 
magnetic braking in post-binary interaction products. Similarly 
magnetic braking may also have contributed to the larger aver¬ 
age spin rate observed in the short- compared to the long-period 
systems. To test this scenario, further higher resolution spectra 
would be necessary as they would enable the search for mag¬ 
netic spin-down signatures in the form of sub-structures in the 
low-velocity part of the spin distributions. 

5.3. The nature of the high-velocity tail and the single-star 
channel toward long-duration GRBs 


|Ramfrez-Agudelo et al. ( |2013| suggested that the high-t; e sin i 
tail identified in the projected spin-rate distribution of the 
presumed-single O-type stars in the VFTS is populated by binary 
interaction products that have been spun up through mass trans¬ 
fer or coalescence. A test of this scenario is to determine the pro¬ 
jected spin-rate distribution of primaries and to check whether a 
high-velocity tail is present or not. If the tail in the sample of 
presumed-single stars is dominated by post-interaction products 
one may anticipate that it is less pronounced or even absent in 
the binary sample. Indeed tides in pre-interaction binarie s can 
not spin the components up beyond 300 km s -1 (see Sect. 5.1 ). 
We find that 2.6+1.8% of the primaries spin faster than 300 


km s -1 , compared to 10.4+1.9 % in the presumed-single sample. 
This supports the hypothesis that the high-velocity tail is actually 
dominated by post-interaction products and that genuine single 
stars spinning faster than 300 km s -1 are rare and maybe absent 
altogether. 

Rapidly spinning single stars have been considered as a po¬ 
tential channel tow ards the production of long-duration gamma- 
ray b ursts (GRBs; |Yoon & Langer||2005| |Woosley & Heger 
2006). If stars are bom with rotational speeds of at least 300- 
400 kms -1 , depending on metallicity, they will evolve quasi- 
chem ically homogeneously (e.g., Maeder 1987; |Brott et ah 
2011 j), i.e., they remain compact throughout their hydrogen 


burning phase and become Wolf-Ray et stars directly afterwards 
In this way their stellar cores may retain enough angular mo¬ 
mentum to fulfill the requirements of the collapsar model to¬ 
wards GRBs ( |Woosley||1993| ), especially in low-metallicity en¬ 
vironments where stellar winds of Wolf-Rayet stars are weaker, 
hence remove less angular momentum ( |Vink & de Koter|2005| . 

If rapidly spinning single stars are rarely formed or not 
formed at all this particular channel toward GRB progenitors 
may not occur in nature. It is however too early to tell whether 
or not this implies that single stars may not produce GRBs al¬ 
together, as stars with low magnetic torques ( ]Woosley|2011|) or 
core collapse early-on in the Wolf-Rayet phase ( [Vink et al.|201 1| 
Grafener et al.|2012| ) might provide alternative routes. Spun-up 
mass gainers in close binaries might evolve chemically homo¬ 
geneously and provide a GRB at the end of their fife ( Cantiello 
et al.|2007| , such that the total number of GRBs predicted from 
chemically homogeneous evolution does not necessarily change. 


5.4. The nature of the low-velocity peak and the formation of 
massive stars 

The low-ARV sources are mostly SB1 systems with relatively 
long orbital periods (with probable periods -10-1000 d) and/or 
low mass ratio systems, neglecting inclination and sampling ef¬ 
fects. They are likely dominated by pre-interaction systems that 
do not undergo significant tidal synchronization (Fig. ED- The 
few systems that are located in the region of the parameters 
where tidal interaction is expected to take place (shaded ar¬ 
eas in Fig. ED occupy the lower left corner, where the tidal 
synchronization timescale is the longest. Moreover, the stellar 
winds of the bulk of this sample are too weak to induce signif¬ 
icant angular-momentum loss, and spin-down due to envelope 
expansion as a result of secular evolution may be neglected. This 
points to the idea that the spin rates of the low-ARV, wide and 
low-q spectroscopic binaries have seldom been modified since 
the stars were formed and thus reflect the outcome of the forma¬ 
tion process. 

While the formation process of massive stars is heavily de¬ 
bated (e.g., |Zinnecker & Yorke|2007||Tan et ah 2014[|Krumholz 
2014 ), most theories agree on the need for disk-mediated ac¬ 
cretion. In its simplest form, the collapse of the initial natal 
cloud and the gas-accretion phase concentrate by far more an¬ 
gular momentum in the central region than can be stored in a 
single star. Therefore, the spin rate of genuinely single massive 
stars may naively be expected to be close to the star break up ve¬ 
locity (the so-called angular momentum problem). Gravitational 
torques between star and disk appear capable of limiting the stel¬ 
lar _spin rates to about half the break-up speed v cr \ t (Lin et al. 
2011). This value is however still substantially higher than the 
0.1-0.3 v CT it that is typical for the observed low-velocity peak 


(e.g.,|Conti & Ebbets|1977[|Huang & Gies|2006 JPenny & Gies 


2009]" Ramirez-Agudelo et al.||2013 ). Magnetic breaking might 
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aid in a further spin-down of the stars, either by self-induced 
fields or by fossil fields. Alternatively, the formation of mas¬ 
sive stars as part of a multiple system may provide a reservoir 
to store angular momentum from the natal cloud into the orbit of 
the components. This too may alleviate the need to have a large 
initial stellar spin. 

The low-velocity peak of the presumed-single stars studied 
by [Ramirez-Agudelo et al.| ( 2013] ) may also reflect the initial 
conditions. Here we compare the low-p e sin i regime of the two 
samples to assess whether the spin distribution of presumed- 
single stars is similar to that of low-ARV primaries in binary 
systems. A Kuiper’s test selecting the stars having v e sin i < 170 
km s -1 in both samples fails to identify significant differences 
and we henceforth accept that these samples may be drawn from 
the same parent distribution. 

The resemblance of the low-velocity peak in the v e sin i dis¬ 
tributions of the long-period and/or low-q spectroscopic binaries 
and the presumed-single stars thus suggests similarities in their 
formation process. This could be the case if primordial mag¬ 
netic fields control the birth spin properties of massive stars, be 
they either single or part of a binary system. If the initial spin 
properties are controlled by multiplicity, i.e., by a formation pro¬ 
cess that proceeds through the creation of pairs or higher order 
multiples, the similarity in the spin distributions of presumed- 
single stars and primaries might suggest that true single stars do 
not actually exist. Interestingly, a high angular resolution survey 
of galactic massive stars combined with spectroscopic results 
from the literature has revealed an extreme multiplicity fraction 
(>90%), suggesting that all massive stars may form as part of a 
binary or higher order system and that higher hierarchical sys¬ 
tems are not uncommon (Sana et al. 2014). Therefore, genuine 
single massive stars may not exist. 


the presumably wider systems with lower ARVs. 

3. About one third of our binaries are so close that the spins 
of their components are expected to be tidally synchronized 
with the orbital motion. The current spin rates in such 
close binaries are thus expected to be independent of their 
birth spins. Tidal synchronization may be responsible for 
the larger average spin rates that we observe between 
components of close binaries and that of single stars and 
of stars in longer period systems. A part of this shift might 
also result from magnetic braking affecting a fraction of the 
single-star and long-period binary population. The latter 
scenario would be consistent with magnetic fields being 
generated by strong binary interaction in an earlier epoch of 
evolution of these systems. 

4. The SB2 systems in our sample mostly show radial ve¬ 
locity variation amplitudes above 200 kms" 1 and have 
tidal synchronization timescales that are shorter than their 
main-sequence timescale. On average, we find that the 
Vq sin i values of SB 2 primaries and secondaries are similar, 
as expected for tidally synchronized binaries. 

5. Among the primaries of our binary systems, stars spinning 
faster than p e sin i = 300 km s -1 are about four times as rare 
as in the single star population, in relative terms. Thi s is c on¬ 
sisten t with the hypothesis, proposed b y |Ramirez- Agudelo 
|et al.| ( |20T3] ) and |de Mink et al.| ( |2013| ), that this tail in the 
(presumed) single star distribution is dominated by post¬ 
interaction binary products. In particular, the population of 
wide binaries shows hardly any components spinning faster 
than 300 kms -1 suggesting that massive star formation 
rarely produces such systems, if at all. 


6. Conclusions 


Based on a sample of presumed-single stars in the 30 Doradus 
region in the LMC presented by [Ramirez-Agudelo et al. ( 2013] ) 
we have constructed a calibration that allows to estimate the pro¬ 
jected spin rates of O-type stars from FWHM measurements. We 
have applied this v t sin /-FWHM calibration to determine the 
projected spin rates of the components of 83 SB1 and 31 SB2 
O-type binary systems in 30 Dor. 

Given the uncertainties in individual measurements (20- 
40 kms -1 ) we have focused on the spin-rate distribution, and 
have compared it to that of presumed-single stars. The sample 
is large enough to define sub-samples that could be compared to 
each other. This has lead to the following conclusions: 


1. The overall spin distribution of the O-type binary compo¬ 
nents qualitatively resembles that of the single O stars in 
the 30 Dor region (see |Rarmrez-Agudelo et al.|2013| ). Both 
are dominated by a low-velocity peak (< 200 km s -i ) which 
contains -80% of the samples, and both show a shoulder at 
intermediate velocities (200 km s -1 < v Q sin / < 300 km s -1 ). 
However, the spin distribution of presumed-single stars fea¬ 
tures a clearly identifiable high-velocity tail (> 300 km s -1 ), 
which is hardly populated in the binary population. 


6. The low-velocity peak in the sample of primaries in wide 
spectroscopic binaries, that orbit their companion at dis¬ 
tances of the order of 1-10 AU, is very similar to that of 
the presumed-single stars. This, together with the overall 
qualitative similarity in the spin distribution of single stars 
and binary components, and plausible scenarios for even ex¬ 
plaining the quantitative differences through binary evolu¬ 
tion, suggests that the spins of massive stars are set indepen¬ 
dently of whether they form in single stars or binary systems 
and are controlled by similar physics. 


The results presented in this work focus on the overall spin 
distribution of our binary O-star sample. A detailed analysis of 
the properties of multiple systems - yielding periods rather than 
A RV as proxies of these periods - and their relation to stellar pa¬ 
rameters is a logical next step in assessing the role of tides and 
mass-transfer in binary systems and the differences and similari¬ 
ties in the intrinsic rotational properties of single and binary sys¬ 
tems. Higher spectral resolution data and an emphasis on metal¬ 
lic lines would also offer the possibility to obtain more accurate 
measurements in the low spin regime. These are needed to inves¬ 
tigate possible sub-structures in the low-velocity peak and their 
relation to, e.g., magnetic spin-down. 


2. The main peak of the binary distribution, centered at -100 
kms -1 , is broader and skewed to somewhat higher spin 
rates than in the single star distribution. Similarly, the spin 
distribution of primaries in systems with ARV > 200 km s -1 , 
which we associate with close binaries with P OT \j <10 days, 
is skewed towards higher spin rates than that of primaries in 
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Appendix A: A Bayesian derivation of the intrinsic 
rotation distribution 

A. 1. Analytical form of the v e distribution 

As described in the main text, we adopt a two-component gen¬ 
eral analytical form for the probability density function (PDF) of 
the intrinsic rotational velocity distribution. Following Ramirez- 
Agudelo e t~aL| ( |20 n| , the adopted PDF is composed of a gamma 
distribution and a normal distribution (see Eq.[2]). The five free 
parameters in Eq. [2] (or,/?,//, cr, I y ) also allow us to represent 
PDFs with a variety of shapes, e.g., a velocity peak with varying 
skewness or peak location, and a high-velocity tail that can vary 
from negligible to strong. For example, in the deconvolved rota¬ 
tional velocity distribution of the sample of s ingle O-type stars 
from the VFTS ( [Ramirez-Agudelo et al.|2013] ), the gamma dis¬ 
tribution allowed to represent a low-velocity peak and the nor¬ 
mal distribution was used to model an additional high-velocity 
contribution. 


A. 2. sin i distribution 

We use the intrinsic orbital parameters distributions of the VFTS 
binary sample and the VFTS binary detection probability as a 
function of orbital period, mass-ratio, eccentricit y, samp ling and 
accuracy of the RV measurements obtained in Sana et al. ( [2013 ) 
to compute the intrinsic distribution of the orbital inclinations of 
the detected binaries in the VFTS sample. The obtained inclina¬ 
tion distribution shows a quasi-absence (< 5%) of systems with 
i < 20° and an over abundance of systems with i > 50° compared 
to a distribution computed with r andom orientation of the binary 
plane in the 3D space (Fig. |A.l| ). While the difference between 
the random-orientation distribution and the detected-binary dis¬ 
tribution are real, the overall effect is small enough that we can 
show that - within the quality of our data - our results do not de¬ 
pend of which distribution is adopted. For consistency, we never¬ 
theless proceed by adopting the detected-binary distribution for 
the binary sample and the random-orientation one for the single 
stars. 


A. 3. v Q sin i distribution 

Given a PDF for the intrinsic rotational velocity distribution, 
the associated projected rotational velocity distribution will 
depend on the distribution of orientations of the rotation axis 
with respect to the line of sight. In what follows, we consider 
two different PDFs for the sin i distribution. The PDF for a case 
where the rotation axes are randomly oriented can be written as: 


h( sin i ) = 


1 


(1/ sin 2 i - l) 1 / 2 


(A.l) 


For the sample of primaries of O-type binaries from VFTS, we 
also consider the PDF of sin i computed for th e dete cted spec¬ 
troscopic binary systems in the VFTS (see Sect.|A.2|). Both sin i 


13 












O.H. Ramirez-Agudelo et al.: VFTS XXI: Stellar spin rates of O-type spectroscopic binaries 



sin i 



(Projected) rotational velocity [km s 1 ] 


Fig.A.l. sin z distribution for the two different cases considered: 
a random orientation for the rotation axis, or a rotation axis 
aligned with the binary axis which is computed for the detected 
spectroscopic binaries in the VFTS sample. 


distributions considered are actually very similar, as shown in 


Fig A.l 


The PDF of the projected rotational velocity, q(v e sin z), is 
obtained by convolving the intrinsic rotational velocity PDF with 
the sin z distribution in the following way: 

q(v e sini) = f f(v') / ; / f ’ c Sln 7 j — dv' (A.2) 

\ V' ) V' 

where / is as defined in Eq.[2] and h can be either of the two sin z 
PDFs discussed above. The likelihood of an individual v Q sin z 
measurement given the above model and a measurement uncer¬ 
tainty <r v is then obtained by the convolution of the above PDF 
and a Gaussian distribution (assuming normally distributed er¬ 
rors): 



4>bs(fe shl /) 


»-r 


1 


V27T < 


■ exp 


—(v e sin if 
2cr 2 , 


q(v e sin i-V)dv', (A 


(Projected) rotational velocity [km s 


Fig. A.2. PDF of the intrinsic (red), projected (green), and ob¬ 
servable (blue) rotational velocity distribution for a model with 

a = 4.82, p = 1/25, jj = 205 km s _1 , cr 2 = (l 90 km s _1 ) , I y = 

0.43, and assuming cr v = 20 km s -1 . Upper panel: sin z distribu¬ 
tion assuming random orientation of the rotation axes. Fower 
panel: sin z distribution computed for the detected spectroscopic 
binaries in VFTS is assumed. 


where q is as defined above. 

To illustrate the effect of inclination and measurement un¬ 
certainties on the rotational velocity distribution, we show in 
Fig. |A.2| the intrinsic, projected, and observable velocity dis- 
tributi on for a chosen set of parameters that were identified by 
|Ramfrez-Agudelo et al. ( 2013| as providing the best representa¬ 
tion of the rotational velocity distribution of the single O-type 
stars in VFTS (a = 4.82, ft - 1/25, p - 205 kms" 1 , cr = 
190 km s -1 , I y - 0.43). For this illustration, we assume a mea¬ 
surement uncertainty of <x v = 20 km s -1 . We show an example 
adopting each of the sin z distributions discussed above. As ex¬ 
pected, in both cases, the projected rotational velocity distribu¬ 
tion is shifted to lower velocities due to the effect of inclination, 
and observational errors slightly broaden the distribution. 


In a Bayesian framework, the posterior probability density can 
be written as: 

pm» = \ p(P\®)pm, ( a . 4 ) 

where p(D\ 0) is the likelihood function and /?(0) is the prior 
distribution. The normalization Z = p(D) is independent of 0 
for a given choice of the form of the generative model, so for the 
problem that we are interested in here, we can simply sample 
from p(®\D) without computing Z. Our likelihood function is 
simply the product of individual likelihoods as given by Eq. |A.3[ 
i.e., 


A.4. Bayesian analysis and posterior sampling 

We are interested in the posterior probability density function of 
the model parameters (0 = {a, /3 , p, cr, I y J) given the data (D: 
a set of v e sin z measurements with associated uncertainties cr v ). 


p(D\&) = [ / ob s,,(4e sin ij, cr „ a, /?, p, cr, l y ), (A.5) 

i 

where v Q sin i t corresponds to each individual measurement with 
uncertainty cr vi . We assume uniform priors for all five model 
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a /? H a / 7 


Fig.A.3. One- and two-dimensional projections of the posterior probability distributions of the rotational velocity distribution 
parameters for the sample of primaries of O-type binaries, assuming the sin i distribution computed for the detected binaries in 
VFTS. 


parameters over the following ranges: 0.1 < a < 10, 0 < j3 < 
0.4, 0 < ji < 500, 30 < cr < 450, 0 < I y < 1. 

Finally, to sample efficiently from the posterior distribution 
and obtain a sampling approximation to the posterior PDF, we 
use the Python implementation of the affine-invariant ensem¬ 
ble sampler for Markov chain Monte Carlo (MCMC) available 
through the emcee code (Foreman-Mackey et al. 2013| ). We use 
100 walkers and run the chains for 3500 steps with a burn-in 
phase of 500 steps. Our final results for the parameters of the 
rotational probability density functions for the primaries and for 
the single stars are given in Table [4] and illustrated in Fig. [6] 
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Fig.A.4. One- and two-dimensional projections of the posterior probability distributions of the rotational velocity distribution 
parameters for the sample of apparently single O-type stars, assuming a random orientation of the rotation axis. 
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